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ABSTRACT 


y 

We  study  the  optimal  control  of  a shuttle  system  consisting  of  a 

single  infinite  capacity  carrier  transporting  passengers  between  two 

/ i 

terminals.  Passengers  arrive  according  to  independent  Poisson  processes, 
and  dispatching  decisions  to  hold  the  carrier  for  more  passengers  can  be 
made  at  only  one  of  the  terminals.  The  objective  is  minimization  of 
the  long-run  average  of  a linear  passenger  waiting  cost  and  a fixed 
charge  per  trip  made.  When  complete  information  about  the  system  state 
is  available,  and  travel  times  are  deterministic,  we  prove  that-^it  is 
optimal  to  dispatch  the  carrier  if,  and  only  if,  the  total  number  of 
passengers  waiting  at  both  terminals  is  greater  than  a cutoff  value. 

An  iterative  method  for  computation  of  the  cutoff  value  is  given  and 

i 

-w«a— finest  hat  it  can  be  well  approximated  by  a function  of  system  costs 
and  parameters  similar  to  the  economic  lot  size  formula.  A (possibly 
non-optimal)  dispatching  rule  is  proposed  for  the  case  when  only  the 
number  of  passengers  waiting  at  one  terminal  is  known,  and  its  efficiency 
is  compared  to  that  of  the  aforementioned  optimal  rule.  Extensions  to 
other  optimality  criteria  and  to  the  case  of  stochastic  travel  times 
are  outlined. 
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INTRODUCTION 


In  this  paper  we  study  a mathematical  model  of  a transportation  system 
in  which  passengers  (or  goods)  are  transported  between  two  terminals  by  a 
single  carrier  with  very  large  capacity  compared  to  the  arrival  rate  of  pas- 
sengers. Our  goal  is  to  determine  optimal  operating  rules  for  such  systems 
In  which  dispatching  decisions  are  made  at  only  one  of  the  terminals. 

The  system  we  study  has  the  following  characteristics: 

(1)  There  are  two  terminals  labeled  1 and  2 at  which  passengers  arrive 
individually  according  to  independent  Poisson  processes  with  ar- 
rival rates  and  A , respectively.  All  arriving  passengers  wait 
to  be  transported  to  the  other  terminal  where  they  exit  the  system. 

(2)  There  is  a single  carrier  which  shuttles  back  and  forth  between 
the  terminals.  The  carrier  has  an  infinite  capacity,  meaning  that 
on  departure  from  the  terminal  it  can  accommodate  all  the  waiting 
passengers.  For  the  majority  of  our  analysis,  we  assume  that  inter- 
terminal travel  time  is  deterministic  and  independent  of  the  load 
carried.  Some  analysis  when  travel  times  are  random  variables 
will  be  found  in  Section  IV. 


Any  views  expressed  in  this  paper  are  those  of  the  authors.  They 
should  not  be  interpreted  as  reflecting  the  views  of  The  New  York  City- 
Rand  Institute  or  the  official  opinion  or  policy  of  the  City  of  New  York 
Papers  are  reproduced  by  The  Rand  Corporation  as  a courtesy  to  members 
of  its  staff. 


-2- 


(3)  The  shuttle  is  operated  with  the  following  dispatching  mechanism : 

At  one  of  the  terminals,  the  carrier  never  waits.  Rather,  it 
visits  this  station,  instantaneously  picking  up  those  passengers 
who  are  already  there — much  as  a bus  does  at  the  typical  urban  bus 
stop.  At  the  other  terminal,  the  carrier  can  wait.  In  fact,  the 
only  control  possible  is  exerted  by  deciding  how  long  it  should 
wait  there.  Our  concern  in  this  paper  is  with  the  determination 
of  optimal  rules  for  dispatching  the  carrier  at  the  terminal  where 
it  can  wait. 

We  use  minimum  long-run  average  cost  per  unit  time  as  our  standard  of 
optimality  and  consider  two  kinds  of  costs:  a fixed  cost  of  making  a trip,  K, 

and  a cost  of  passenger  waiting,  h,  per  passenger  per  unit  time.  Our  principa' 
result,  developed  in  Section  I,  is  that,  when  the  dispatcher  has  available 
complete  information  about  the  state  of  the  system,  there  is  an  optimal  dis- 
patching rule  of  the  form:  The  carrier  should  depart  if, and  only  if, the 
total  number  of  customers  waiting  in  the  system  exceeds  a constant.  Our 
main  tool  is  a result  of  Deb  and  Serfozo  [2]  for  a one  terminal  system. 

We  refer  to  such  a rule  as  a control  limit  rule  and  to  the  constant 
as  a control  limit.  An  expression  for  average  cost  as  a function  of  the 
control  limit  constant  is  given  in  Section  II.  Computations  with  this 
expression  indicate  that  the  optimal  value  of  the  control  limit,  k* , is  well 
approximated  by  a simple  function  of  the  cost  and  arrival  parameters, 
namely  ' 

k*  = /2(X1  + X2)K/h  . 

Our  computations  also  indicate  that,  if  one  can  choose  which  terminal  to 
dispatch  from,  average  cost  is  slightly  lower  if  the  terminal  with  more 
passenger  arrivals  is  chosen. 

In  Section  III  we  consider  the  case  where  the  dispatcher  does  not  know 
how  many  passengers  are  waiting  at  the  other  terminal.  The  results  of 
Section  I suggest  the  desirability  of  a rule  of  the  form:  The  carrier 
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should  depart  from  the  terminal  if,  and  only  if,  the  number  of  customers 
waiting  at  the  terminal  plus  a multiple  of  the  time  the  carrier  has  waited 
there  exceeds  a control  limit.  An  expression  for  average  cost  as  a function 
of  the  multiple,  call  it  3,  and  the  control  limit  is  obtained,  and  calcula- 
tions reveal  that  the  optimal  value  of  3 is  typically  larger  than  the  pas- 
senger arrival  rate  at  the  other  terminal.  More  important  is  the  observa- 
tion that,  for  a fixed  trip  rate,  the  dispatching  rule  using  only  passenger 
waiting  information  from  one  terminal  results  in  25  percent  to  75  percent 
higher  average  queue  sizes  than  the  optimal  dispatching  rule  using  complete 
passenger  waiting  information. 

Section  IV  is  devoted  to  stochastic  travel  times  where  it  is  indicated 
that  the  result  of  Section  1 still  holds.  In  Section  V,  an  alternative  way 
of  viewing  the  problem — minimizing  the  average  number  of  passengers  waiting 
at  the  terminals  subject  to  a constraint  on  the  average  number  of  trips  per 
hour  by  the  shuttle — is  considered.  It  is  shown  that,  for  properly  chosen 
values  of  the  trip  rate  constraint,  rules  of  the  form  shown  optimal  for  the 
problem  of  Section  I are  again  optimal. 

One  may  question  the  relevance  of  these  results  since  real  shuttles 
have  finite  capacity  and  it  would  seem  that,  for  such  systems,  one  ought 
then  to  pay  attention  t^  i^  and  i^  (the  number  waiting  at  terminals  1 and  2, 

respectively)  individually  and  not  just  to  their  sum.  However,  the  infinite 

capacity  model  would  be  useful  if  the  following  conjectures  were  true.  (We 
think  they  are.)  First:  An  optimal  infinite  capacity  rule  is  an  upper 

hound  for  an  optimal  finite  capacity  rule.  That  is,  if  an  optimal  rule  in 
the  infinite  capacity  case  is  to  dispatch  if,  and  only  if,  i^  + i^  k*, 
then  an  optimal  rule  in  the  finite  case  is  to  dispatch  if  i^  + i^  >_  k*, 

and  possibly  in  other  situations  also.  Second:  If  k*  is  less  than  the 

capacity  of  the  shuttle  and  if,  in  the  infinite  case,  the  probability  that 
the  shuttle  has  room  for  all  who  are  waiting  when  it  arrives  at  a terminal 
is  close  to  one,  then  an  optimal  rule  for  the  finite  case  is  close  to: 
dispatch  if,  and  only  if,  i^  + 1^  >_  k*. 

There  has  been  relatively  little  reported  work  on  shuttles  or  on 
other  multiterminal  systems  \ ' th  probabilistic  passenger  arrivals.  Our 
previous  work  [3]  was  restricted  to  the  case  of  a shuttle  which  could 
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transport  only  one  passenger  at  a time.  Barnett  [1]  has  compared  average 
number  waiting  for  several  operating  rules  for  an  infinite  capacity  shuttle. 
His  approach  differs  from  ours  in  several  respects:  Ours  has  a cost  per 

trip  and  allows  travel  time  to  be  stochastic,  while  his  allows  the  shuttle 
to  wait  at  either  of  the  terminals. 


I.  THE  FORM  OF  AN  OPTIMAL  DISPATCHINO  POLICY 


For  definiteness,  suppose  that  the  shuttle  may  wait  only  at  terminal  1. 

Let  X^  and  X^  be  the  passenger  arrival  rates,  and  r the  round  trip  travel 

time  (r/2  in  each  direction). 

Decisions  are  permitted  upon  the  arrival  of  the  shuttle  at  terminal  1 
and  upon  every  passenger  arrival — at  either  terminal — thereafter,  as  long  as 
the  shuttle  is  waiting  at  terminal  1.  The  cost  of  a round  trip  is  K > 0, 
Passenger  waiting  cost  is  linear  with  rate  h > 0 so  that,  when  there  are  a 
total  of  i passengers  waiting  at  the  terminals,  cost  is  incurred  at  rate 
h • i. 

The  state  at  epoch  t is  given  by  the  pair  _i(t)  = ( i^ ( t ) , i^Ct)), 

where  t ) is  the  number  of  passengers  waiting  at  terminal  k at  epoch  t. 

Wherever  clarity  will  not  be  lost,  we  will  suppress  the  t and  write 
ij  , i-2’  etc-  Let  * = + *2’  6l  = > and  = ^ ‘ 

Problem  1 

We  will  restrict  ourselves  to  stationary  Markovian  (SM)  control 
policies:  that  is,  those  where  the  decision  made  when  the  system  is  in 

state  i at  a review  epoch  depends  only  on  the  state  at  that  epoch.  A 
SM  rule  can  be  viewed  as  a specification,  for  each  (i^,  ijJt  t*le 
probability  that  we  dispatch  the  shuttle  when  (i  , i0)  are  waiting. 
Consequently,  if  we  let  I = {0,  1,  2,  ...},  then  any  function  A(.,.) 
from  I x I onto  [0,  1]  is  a SM  rule  for  our  system,  and  any  SM  rule 
induces  such  a function. 

Under  any  sucli  operating  rule,  our  system,  when  observed  at  decsion 
epochs,  forms  a semi-Markov  process  whose  key  elements  are  a set  of  tran- 
sition probabilities,  sojourn  times  and  one  step  costs.  We  follow  the 
notation  of  Ross  [5,  Chapter  7].  There  are  2 actions:  the  carrier  can  be 
held  until  the  next  decision  epoch,  denoted  by  W (for  wait);  or  the  carrier 
can  be  dispatched,  denoted  by  G (for  go).  The  transition  probabilities  are: 
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»I,J>6k  ' V*  k ‘ ^ 2 


P.  m = am  ^Ir^am  <Sr/2) 
j^m  m.  1 m„  l 


n^,  m?  = 0, 


“ z i 

where  a.(z)  = a z / j ! is  the  Poisson  probability 


The  sojourn  or  transition  times  are  exponential 
expectation 


(1.1) 

1,  ...  (1.2) 
function . 

random  variables  with 


tW  = 1/X,  (1.3) 

when  the  shuttle  waits.  If  the  shuttle  Roes,  the  transition  time  is  a 
constant  and 


G 

c = r . 


(1.4) 


We  shall  be  interested  in  long-run  average  cost  per  unit  time  and  so 
must  concern  ourselves  with  the  expected  transition  costs; 

If  the  shuttle  waits,  expected  transition  cost  is  given  by 

CW(i)  = h • (ij  + i2)/X,  (1.5) 

since  the  number  of  passengers  waiting  is  constant  over  the  interval  until 
the  next  one  arrives.  To  obtain  expected  transition  cost  if  the  shuttle 
goes,  we  define  the  random  variable  I^(t|iL)  as  the  number  of  passengers 
waiting  at  terminal  k t minutes  after  the  shuttle  departs  from  terminal  1, 
if  the  state  was  i^  just  before  departure,  (t  <_  r).  When  the  shuttle  leaves, 
all  passengers  now  waiting  at  terminal  1 are  picked  up,  and,  after  r/2,  all 
then  waiting  at  terminal  2 are  picked  up.  Additions  to  both  queues  are  by 
Poisson  arrivals.  So  we  have;  I^(t|i_)  is  a Poisson  random  variable  with 
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parameter  A^t,  and  ^(tljL)  - i,  is  Poisson  with  parameter  X2t  for 
0 ^ t ^ r/2,  while  1 2 ( c I i)  is  Poisson  with  parameter  X2(t  - r/2)  for 
r/2  t <_  r.  Consequently,  we  have 


CG(i)  = K + 


X2  i 2 

Xi  + if1-  + — 
12  r 


(1.6) 


We  wish  to  find  a dispatching  policy  that  minimizes  the  average 
cost  per  unit  time  of  operating  this  svstem. 

It  is  well  known  that  the  average  cost  per  unit  time  of  operating 
under  control  rule  A is  given  by 


f?1(A) 


V tt.  • (A(i)C°(i)  + [1  - A(i)]CW(i)) 


l tt±  • (A(i)rG  + [1  - A(i)]iW) 


(1.7) 


where  is  the  stationary  probability  distribution  of  the  imbedded 

Markov  chain  under  policy  A.  (The  subscript  1 refers  to  Problem  1.) 

If  these  probabilities  do  not  exist  (implying  that  the  chain  is  transient 
or  null  recurrent),  we  take  g^(A)  = •*.  Of  course,  for  any  policy  with 
A(^L)  = 1 for  i^  _>  i^*  and  i 7 _>  i9*,  the  process  is  ergodic  and  g^(A)  will 
be  finite. 

It  is  also  well  known  (Ross  [5])  that  a policy  which  minimizes  g^ 
results  from  choosing  the  minimizing  action  for  each  i_  which  satisfies 
the  functional  equation: 

’ 

„W  ...  W r W 

[c  (i>  - ST  + l vj 

CC(i)  - gT°  + I Pjj  V. 


v^  = min 


(1.8) 
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However,  (1.8)  is  inconvenient  in  several  ways — it  does  not  reveal  the 
form  of  the  optimal  policy  and  computations  cannot  be  made  directly  with 
it  since  there  are  an  infinite  number  of  equations  and  variables. 

In  order  to  show  that  an  optimal  SM  policy  exists  for  Problem  1 of 
the  form:  Dispatch  the  shuttle  if  (and  only  if)  ^ + i9  _>  k,  and  also 

to  determine  a value  of  k which  minimizes  average  cost,  we  consider  two 
related  problems.  The  first  of  these  problems  is  used  to  show  that  the 
optimal  policy  need  depend  only  on  knowledge  of  i^  + i and  does  not 
require  knowing  i^  and  i.,  individually. 


Problem  2 

This  problem  is  the  same  as  Problem  1 except  that  we  also  charge  a 
waiting  cost  for  passengers  from  terminal  1 during  their  trip  to  terminal  2. 
Formally,  we  then  have 


CG(i) 


(1.9) 


and  Problem  2 is  to  minimize  average  cost  for  the  system  given  by  (1.1)- 
(1.5)  and  (1.9). 

Consider  any  SM  policy  A.  It  is  elementary  to  prove  that  the  long-run 
average  cost  obtained  when  this  rule  is  applied  to  system  2 is  related  to 
the  long-run  average  cost  when  it  is  applied  to  system  1 by 


g2(A)  = gl(A)  + h.^r/2,  (1.10) 

because  every  passenger  from  terminal  1 eventually  travels  to  terminal  2. 

Equation  (1.10)  implies:  any  SM  policy  that  is  optimal  for  Problem  1 

is  also  optimal  for  Problem  2,  and  vice  versa. 


More  formally:  over  an  interval  t units  long,  X^t  + e(t)  passengers 

arrive  at  terminal  1 and  travel  to  terminal  2,  incurring  total  "extra" 
waiting  cost  hr/2(Xjt  + e(t)).  The  average  of  this  cost  per  unit  time 
is  then  hr/2(X^  + e(t)/t),  and  by  the  law  of  large  numbers  e(t)/t  -*•  0 as 

t ->  CO  # 


We  now  argue  that  an  optimal  control  rule  need  depend  only  on  i^  + i 9 . 

A formal  argument,  relating  Problem  2 to  Problem  3 (which  will  be  introduced 
shortly),  is  given  in  the  Appendix.  An  informal  one  goes  as  follows: 

(1)  C1,  (i)  can  be  written  in  the  form  K'  + b(i^  + i9)  and  hence  is 

a function  only  of  i.  + i . 

W x ^ 

(2)  C (i_)  = h/A(i  + i9),  again  a function  only  of  i + i . 

y ^ G ^ 

(3)  T and  T and  p,  do  not  depend  on  _i  at  all. 

W — 

(4)  The  p.  are  such  that  m + m„  = i + i + 1;  that  is,  if 

X ) m X Z X Z 

decision  W is  made,  the  system  evolves  to  a new  state  with 

exactly  one  more  customer  waiting  (somewhere)  in  the  system. 

It  has  been  observed  above  that  the  expected  cost  C^(i)  de- 

W G 

pended  only  on  the  sum  i^  + i , and  the  costs  C (m)  or  C (m) , 
one  of  which  is  now  to  be  incurred,  depend  only  on  the  sum 
ij  + i0  + 1 so  that  all  future  evolution  and  costs  of  concern 
to  us  are  dependent  only  on  i^  + i9. 


Problem  3 

Deb  and  Serfozo  [2]  consider  an  infinite  capacity  batch  service  queueing 
system  with  a Poisson  arrival  process  and  a similar  cost  and  decision  structure 
to  that  of  Problem  1.  They  show  that  long-run  average  cost  per  unit  time  is 
minimized  by  a policy  of  beginning  a service  only  if  the  number  of  customers 
in  queue  is  at  least  k*.  Their  problem  is  similar  to  our  Problem  1,  but  their 
system  takes  into  account  only  activity  at  terminal  1.  We  will  show  how  we  can 
relate  their  problem  to  Problem  2,  and  hence  to  Problem  1.  First,  we  describe 
their  problem:  Consider  a single  terminal  (Deb-Serfozo)  system  with  arrival  rat 

A^  + ^2  ~ A at  the  single  terminal.  Suppose  that  round  trip  travel  time  is 
deterministic  and  equal  to  r,  that  B + dj  is  the  cost  of  a round  trip  made 
with  j passengers  on  board,  and  h > 0 is  the  waiting  cost  rate.  Let  i denote 
the  number  of  waiting  passengers.  Decisions  to  wait  or  go  are  permitted  at 
the  epochs  of  the  arrival  of  the  carrier  at  the  terminal  and  at  passenger 
arrival  epochs  as  long  as  the  carrier  is  waiting  at  the  terminal.  This  is 
a semi-Markov  decision  process  with: 
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CC(i)  •=  B + di  + — - X 


(1.11) 


C"(i)  = hi/\ 


(1.12) 


Pi,i+1  = 1 


Pij  = aj(Ar) 


= 1/X 


(1.13) 


(1.14) 


(1.3) 


G 

t = r 


(1.4) 


As  before,  we  wish  to  minimize  the  long-run  average  cost  per  unit 
time  for  this  system,  and  Deb  and  Serfozo  have  proved  that,  among  the 
optimal  rules,  there  is  a rule  such  that  the  carrier  is  dispatched  if, 
and  only  if,  the  number  of  passengers  waiting  in  line  exceeds  a minimum 
level  called,  say,  k*. 

Now,  recall  our  assertion  that  for  Problem  2 we  need 
only  concern  ourselves  with  the  total  number  of  passengers  waiting. 

That  means  that,  for  the  2-terminal  system,  we  need  only  observe  the 
total  number  in  the  system,  say,  i.  Problem  2 when  viewed  in  this  way 
forms  a semi-Markov  process  whose  key  parameters  are  defined  by  (1.3), 
(1.4),  (1.9),  (1.12),  (1.13)  and  by 

j 

pij = l=Q  ak(xir)  aj-k(x2r/2)  = yv  v r)-  (i-i5) 

2 

Letting  d = hr/2  and  B = K.  - X hr  / 4 yields  equality  between  (1.9)  and 
(1.11)  so  the  only  remaining  difference  between  Problems  2 and  3 is  the 
difference  between  (1.14)  and  (1.15).  But  the  particular  form  of  (1.14) 
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is  incidental.  Examination  of  Deb  and  Serfozo’s  proofs  reveals  it  is  suf- 

Q 

ficient  for  the  p to  be  independent  of  i.  Since  this  is  also  the  case 
for  (1.15)  we  claim: 

Theorem  1 

There  exists  an  optimal  (long-run  average  cost)  dispatching  rule  of 
the  form:  Dispatch  the  carrier  if,  and  only  if,  the  total  number  of 

passengers  waiting  in  the  system  exceeds  a control  limit  k*. 

Note  that  this  result  can  be  extended  to  more  than  two  terminals  as 
long  as  the  shuttle  is  held  only  at  terminal  1. 
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II.  FINDING  AN  OPTIMAL  DISPATCHING  NUMBER 

Having  shown  that  an  optimal  SM  rule  is  of  the  form:  Dispatch  if, 
and  only  if,  i^  + i9  _>  k,  we  present  a method  by  whicli  the  optimal  k can 
be  determined  given  the  system  parameters  and  costs.  We  provide  below 
an  expression  for  g,  the  average  cost  per  unit  time, as  a function  of  k 
which  can  easily  be  evaluated  with  the  aid  of  a computer.  We  have  made 
calculations  with  this  expression  showing  how  the  optimal  k depends  on 
the  parameters  (A^,  X9,  r,  h,  K) . These  results  will  be  discussed  at  the 
end  of  this  section. 

We  do  our  calculations  using  Problem  3,  and  then,  using  (1.10), 

subtract  hA^r/2  from  the  resulting  value  of  g to  get  the  average  erst 

for  Problem  1.  Following  Deb  and  Serfozo,  let  T^  be  the  epoch  of  the 

nth  shuttle  arrival  at  terminal  1,  and  Tq  0.  Let  u = (\ ^ + X^/2)r. 

Let  X be  the  total  number  of  passengers  waiting  at  both  terminals  at  epoch 
n 

T^.  {X^}  is  a sequence  of  independent  identically  distributed  random 

variables  with 

P {X^  = j}  = a^ (u)  = e 3u^/j!. 

Let  Z.  be  the  expected  cost  in  the  interval  (T  , T given  X = i. 
j n n+i  n 

Then,  analogous  to  (1.7),  we  have  (writing  g^(k)  to  denote  the  dependence 
of  the  average  cost  for  Problem  3 on  the  control  limit  value  k) 


(2.1) 

I 


g3(k) 


l ^ («)z1 


|.0  »]<“>  • E<Tn+l  - TJXn  ‘ 
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From  (1.3)— (1.4)  we  have: 


E (T 


- T |X  = j)  = 

n 1 n J ' 


r If  j >_  k 
r + (k  - j ) / X 


if  j < k. 


Using  a process  analogous  to  Ik(t|i)  from  Section  I and  omitting  details: 
If  j > k,  Zj  = K + h(|£  + jr/2). 

If  j < k,  Zj  = K + h|j  (k  - j ) /X  + (~  + kr/2)  + (k  - j - 1)  (k  - j)/2X 


Rewriting  the  kr/2  term  as  jr/2  + (k  - j ) r / 2 and  observing  that 

a (a)hjr/2  = , 

j=0  J 


K + h 


R3(k) 


k-1 

ur  + (1/2A)  l (k  - j)0  + k + \r  - l)a  fa) 
i=Q  1 j 


k-1 

r + (1/D  l (k  - j)a  (a) 
j=0  1 


(2.2) 


and 


8r(k)  = g3(k)  - h^r/2. 


(2.3) 
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We  have  made  computations  with  (2.3)  for  various  values  of  A^,  , and  K. 

We  have  set  r = h = 1,  and,  in  so  doing,  effectively  have  chosen  the  time 
scale  so  that  r = 1 and  cost  units  so  that  h = 1.  We  have  computed  g^(k) 
for  k = 0,  1,  2,  ...,  30  for  490  cases:  7 values  of  K,  namely  0,  i,  2,  4, 

8,  16,  and  32;  for  each  value  of  K,  7 values  of  X and  7 of  u(where 
u = (A^  + \^/2)r  = A^  + A.,/2)  , namely  .25,  .5,  1,  2,  4,  8,  and  16; 
and,  for  each  value  of  A and  u,  5 values  of  A^.  For  each  value  of  A, 
the  values  of  A^  are  .125A,  ,25\,  .5A,  .75A,  and  .875A;  similarly, 
for  each  value  of  a,  the  values  of  ^ are  .125u,  ...,  .875u.  In  each 
case,  we  find  that  g(0)  > g(l),  implying  tiiat  it  is  always  a good  idea 
to  wait  for  at  least  one  passenger,  and  that  g(k)  is  convex  in  k on 
{1,  2,  3,  ...»  30).  The  minimizing  value  of  k,  call  it  k*,  increases 
with  A and  with  K,  and,  as  Table  1 indicates,  to  a first  approximation 
k*  is  equal  to  /2AK.  (This  approximation  is  suggested  by  Deb  and 
Serfo2o's  result  for  a variant  of  Problem  3.)  A complete  table  of 
results  appears  in  the  Appendix. 

To  see  the  effect  of  changing  h,  we  note  that  doubling  h would 
double  the  value  of  each  cost  unit  and  thus  halve  the  value  of  K. 

To  see  the  effect  of  changing  r,  we  note  tiiat  doubling  r would  double 
the  size  of  each  time  unit  and  hence  double  A^  and  A ^ ; h,  being  the 
holding  cost  per  customer  per  unit  time,  would  also  double,  thus 
halving  K.  As  a result,  the  rule  of  thumb  that  k*  is  equal  to 
>2 AK  when  r = h = 1 becomes  k*  is  approximately  equal  to  /2AK/h 
wiien  r and  h are  unconstrained. 


III.  A RULE  WHEN  THE  QUEUE  SIZE  AT  TERMINAL  2 IS  NOT  KNOWN 

So,  when  the  number  of  passengers  waiting  at  each  terminal  is  known  to 
the  controller,  it  is  optimal  to  dispatch  the  shuttle  if,  and  only  if,  the 
total  number  of  passengers  waiting  at  both  terminals  is  at  least  k.  We  have 
not  been  able  to  derive  the  form  of  an  optimal  dispatch  policy  when  the  con- 
troller has  knowledge  only  of  the  number  of  passengers  waiting  at  terminal  1. 
Instead,  we  suggest  a family  of  rules  which  are  likely  to  be  quite  good. 

These  rules  are  motivated  by  the  form  of  the  rule  when  complete  passenger 
waiting  information  is  available  — that  is,  go  if  i^  + _>  k.  In  the  situa- 

tion at  hand,  i9  is  not  known,  but,  if  the  carrier  has  been  waiting  at  ter- 
minal 1 for  t minutes,  he  left  terminal  2 empty  t + r/2  minutes  ago  anil, 
hence,  the  expected  number  of  passengers  waiting  at  terminal  2 is  A^(t  + r/2), 
and  we  consider  dispatch  rules  of  the  form:  Go  if  i^  + B(t  + r/2)  k. 

Dropping  the  subscript,  noting  that  k need  not  be  integer,  and  defining 
u = k - 6r/2,  the  rule  can  be  expressed  as:  Go  if  i + tit  ^ ii.  Equations 

(3.5)-(3.9),  which  are  developed  below,  permit  calculation  of  average  cost 
per  unit  time,  trip  rate,  and  average  number  of  passengers  waiting  as  func- 
tions of  the  pair  (u,  8) . Calculations  made  with  these  equations  indicate 
that  the  optimal  pair  (u*,  d*)  is  typically  such  that  C*  is  considerablv 
larger  than  A^.  For  example,  for  A^  = = .5  and  r = h = K = 1, 

(u*,  ti*)  = (.7,  .66),  so  d*  is  about  30  percent  larger  than  A . 

More  important,  these  calculations  coupled  with  those  described  in  the 
previous  section,  suggest  the  value  of  knowing  the  number  waiting  at  ter- 
minal 2.  In  this  discussion,  it  seems  more  informative  to  make  comparisons 
of  the  average  number  waiting  rather  than  average  cost.  If  an  optimal 
SM  rule  for  Problem  1 is  used  as  a point  of  comparison,  we  find  that  the 
best  of  the  (u,  d)  rules  with  approximately  the  same  trip  rate  typically 
has  the  average  number  of  passengers  waiting  25  to  75  percent  higher. 

For  example,  consider  the  parameters  in  the  previous  paragraph.  The 
opt imal  SM  rule  for  this  case  when  both  queues  sizes  are  known  to  the 
controller  is  k*  = 1 and  it  yields  .68  trips  per  unit  time  and  an 
average  of  .42  passengers  waiting.  The  (.7,  .66)  rule  which  is  best 


N'uraber  of  passengers  waiting  at  terminal 


in  the  (u*,  b*)  family  yields  .67  trips  per  unit  time  and  an  average  of 
.62  passengers  waiting.  The  main  effect  of  the  lack  of  information  is 
a 45  percent  increase  in  the  average  number  of  passengers  waiting. 

To  carry  out  our  analysis,  we  again  work  with  the  {T  } and  {X^} 
processes  defined  in  Section  II,  except  that  is  now  the  number  of 
passengers  waiting  at  terminal  1 when  the  shuttle  arrives  there.  Let 
be  the  time  the  shuttle  waits  at  terminal  1 when  X^  = j . Clearly, 

D,  = 0 for  j > u.  For  j < u 


Let  us  begin  by  relating  to  passenger  arrivals  at  terminal  1.  It  is 
helpful  to  start  with  Fig.  1: 


Fig.  1-  Two  realizations  of  the  (u,  b)  dispatching  mechanism 


* 
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If  arrivals  to  terminal  1 occur  at  t^  and  t^  after  the  shuttle  arrival  at 
T , as  illustrated  by  the  solid  line  in  Fig.  1,  then  = t0,  and  we  will 
say  that  the  shuttle  departure  is  "triggered  by  an  arrival."  If  arrivals 
to  terminal  1 occur  at  u^,  u^,  and  u^>  as  illustrated  by  the  dotted  line  in 
Fig.  1,  then  = (u  - 2)/|3,  and  we  will  say  that  the  shuttle  departure  is 
"triggered  by  time." 

Let  A^  be  the  time  until  i arrivals  occur  at  terminal  1;  has  a 

gamma  distribution  with  parameters  i and  X.  From  Fig.  1 we  can  see  that, 

for  time  triggering  to  occur,  there  will  be  exactly  i arrivals  at  terminal  1 

for  i = 0,  1,  ...,  m where  m is  the  greatest  integer  less  than  u and 

the  event  {Dn  = U 0 1 1 = the  event  {exactly  i arrivals  at  terminal  1 in 
. u p 

[0,  ^)). 

For  arrival  triggering,  we  see  that 


(A.  = t} 

l 

for 

t£ 

! u - 

l e 

JL 

u 

’ 6, 

{A.,  = t> 

for 

t£ 

r 

1 u “ 

l 6 

2 

u_ 

9 

r.  .1  ^ _ ftu-m-1]  u-m 

Am+1  = t}  f°r  tE  ( I * 6 


For  Xn=j<u,  it  should  be  clear  that  shifting  the  "dispatching  line"  in  Fig.  1 
down  an  amount  j allows  the  previous  analysis  to  be  repeated.  It  is  sufficient 
then  to  let  u - j play  the  role  of  u,  and  we  obtain  for 


(D, 


4^)  = 


exactly  i arrivals  at 


terminal  1 in  [0, 


u ~ 1 - J> 

e J 


for  i = 0,  1 m-j; 

j = 0,  .... 


m 


(3.1) 
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and 


{D  = t}  = {A±  = t}  for 


te 


riu  - 1 - ji+ 

3 


u - i - i + 1 

3 


for  i = 1,  . . . , m - j + 1; 
j = 0,  . . . , m. 


From  (3.1)  we  get 


P 


D. 

J 


for  i = 0,  . . . , m - j ; 
j = 0,  . . . , m. 


From  (3.2)  we  get 


ti 


> 1 I1"1 

1 

(1  - D! 


— f(Alt)  * 


for 


i J]' 


< t 


u - 1 


* 

Let  C.  be  the  expected  cost  incurred  In  the  interval  (T  , T 


X = j . Then 
n J 


n+1  ■ 


J=0 


*i(iir)cj7? 


j=0 


(3.2) 


(3.3) 


+ 1 


(3.4) 


given 


(A^)  [ED^  + r] . 


(3.5) 
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For  the  denominator,  we  have  ED^  = 0 for  j > m,  and  for  j < m, 


ed  = y 


m-j 


J-'—1  . , 


i=0 


, u - i - i 
A1 8 1 


m-j+1 

+ I 

i=l 


u-i-j+1 

8 

fu-i-il4-  C ‘ Vi-l(Xlt)dt 


i=l 


i 

y 

c 

1 

M. 

1 

u - i - j + l) 

m-j+1  , 

= 1 f 1 

L 

k=0 

ak 

11  B j 

" 3k 

A , ; — ** 

l 1 8 

(3.6) 


+ a . 


X • IL-i-rJ  + .1 
i e 


where,  in  order  to  evaluate  the  integral,  we  note  that  the  integrand  can  be 
rewritten  as  a^(A^t)  • i and  use  the  identity 


fW 

J a (Ax)dx 
1 v n 


T l [a.  (Av)  - a,  (Aw) ] . 
k=0  K K 


(3.7) 


For  the  numerator,  we  have,  omitting  the  details,  for  j > m 


C*  = K + hAr2/2. 


(3.8) 


For  j < in, 


* 

Cj  = A + B + C 


(3.9) 
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A gives  expected  cost  incurred  after  the  shuttle  leaves  terminal  1,  excludin’ 
costs  associated  with  arrivals  occurring  during  the  shuttle's  wait.  B gives 
expected  cost  associated  with  the  shuttle's  wait,  given  time  triggering. 

The  first  of  the  four  terms  multiplied  by  a^  gives  cost  incurred  from 
(1)  those  passengers  waiting  at  terminal  2 when  the  shuttle  arrives  at 
terminal  1,  incurred  during  the  shuttle's  stay  at  terminal  1,  and  (2)  those 
passengers  arriving  at  terminal  2 during  the  shuttle's  stay  at  terminal  1, 
incurred  during  the  shuttle's  trip  to  terminal  2.  The  second  term  gives 
cost  from  arrivals  to  terminal  2 during  the  stay,  incurred  during  the  stay. 
The  third  term  gives  cost  from  those  waiting  at  terminal  1 when  the  shuttle 
arrives  there,  incurred  during  the  stay.  The  fourth  term  gives  cost  from 
arrivals  at  terminal  1 during  the  stay,  incurred  during  the  stay.  C gives 
expected  cost  associated  with  the  shuttle's  stay,  given  arrival  triggering 
In  writing  down  A,  B,  and  C we  have  used  a process  like  I (t|i)  from 
Section  1,  and  have  used  (3.7)  to  simplify  C. 


_ 


For  random  variable  travel  times,  we  have  results  analogous  to  the 
constant  travel  time  case  when  both  queue  sizes  are  known.  This  extension 
is  outlined  briefly  below.  When  dispatching  is  done  without  knowledge  of 
the  situation  at  terminal  2,  a broader  class  of  rules  than  those  investigated 
in  the  previous  section  seems  appropriate  when  travel  times  are  stochastic. 
For  example:  suppose  the  previous  round  trip  travel  time  becomes  known  when 

the  shuttle  arrives  at  terminal  1.  The  longer  it  is,  the  larger  the  number 
of  passengers  likely  to  be  waiting  now  at  terminal  2,  and  it  would  be  natural 
to  include  this  conditional  expectation  in  a rule  for  deciding  when  to  go. 

Let  D and  B be  travel  times  from  terminal  1 to  terminal  2 and  from 
2 to  1 , respectively,  and  let  R = D + B.  We  assume  that  successive  values 
of  I)  are  independent,  with  distribution  function  F,  that  successive  values 
of  B are  independent,  with  distribution  function  G,  and  that  the  D's  and 
B's  are  independent  of  each  other  and  of  the  arrival  and  dispatching  pro- 
cesses. Then,  the  development  parallels  that  of  Section  I,  with  the  fol- 

«» 

lowing  modifications. 

For  Problem  1,  CW,  p'v  and  t'^  are  unchanged  and  C^,  p^  and  are 
given  below. 

We  have  t ' = ER.  (4.1) 

To  get  C*’,  we  look  at  the  random  variable  I^(t|_i)  defined  in  Section  1. 
Let  W^(t)  be  the  expected  waiting  cost  at  terminal  1 conditioned  on  R = t. 
Then , 

Wx(t)  = h fQ  F.I1(y|j.)dy 
= hX1t2/2. 
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Pij 


= l 

jl+J2=1 


where  the  p1.'  are  given  in  (4.3)  does  not  alter  their  validity. 

It  can  be  shown  that  this  replacement  converts  their  problem  into  one  that 
is  equivalent  (in  the  sense  described  in  Section  1 and  in  the  Appendix)  to 
Problem  2 with  stochastic  travel  times.  Given  the  equivalence  of  Problems  1 
and  2 in  this  stochastic  travel  time  case,  we  can  then  conclude  that  a control 
limit  rule  is  optimal  for  Problem  1. 
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Letting  H be  the  distribution  function  of  R,  the  unconditional  expected 
waiting  cost  at  terminal  1 is 


“i  - L <wit: 


/2)dH( t) 


hX^ER  /2 . 


Using  similar  analysis  for  terminal  2,  we  obtain 


C°(i)  = K + hi2ED  + | (X^ER2  + X^.B2)  . 


(4.2) 


Since  arrivals  in  D and  in  B are  independent,  we  can  write 


'ij  = [=0  L aj1-k(A1t)dD(t>  • L ak(xit}  aj2(x2t)dB(t) 


(4.3) 


For  Problem  2,  we  can  add  hijED  to  C (i) , and,  analogously  to  Section  i 
show  that  g9  = g^  + hljED. 

For  Problem  3,  the  transformations  given  in  Section  I are  appropriate 

_W  W , W G ..  n.  _ _G 

tor  C , p and  t . For  x , use  (4.1).  lor  G , use 


CG(i) 


K + (h/2) (X^ER2  + X2EB2)j  + i • hED. 


(4.4) 


And  for  p ’ use 


-ij  = J0  aj (Xt)dH(t) 


(4.5) 


Deb  and  Serfozo  proved  that  a rule  of  the  form  dispatch  the  carrier  if, 
and  only  if,  i > k*  is  optimal  for  the  problem  specified  by  (1.3),  (1.12), 
(1.13),  (4.1),  (4.4),  and  (4.5).  Examination  of  their  proofs  reveals  that 
replacing  (4.5)  by 


d 
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V.  AM  ALTERNATIVE  CRITERION 

An  alternative  to  our  approach  of  assigning  costs  both  to  passenger 
waiting  and  to  trips  made  by  the  carrier  is  to  view  the  minimization  of  pas- 
senger waiting  as  our  objective  while  constraining  the  average  number  of 
trips  made  by  the  carrier  to  not  exceed  a preassigned  limit.  For  any 
rule  A,  let  T(A)  be  the  expected  cycle  time — that  is,  the  expected  time 
between  successive  arrivals  of  the  carrier  at  terminal  1.  Let  W(A)  be 
the  expected  total  waiting  cost  in  a cycle.  Define  0(A)  as  W(A)/T(A). 
Formally,  then,  our  alternative  to  Problem  1 is  to  choose  that  SM  rule  A' 
that  minimizes  Q(A)  over  the  set  of  all  SM  rules  that  satisfy  T(A)  >_  7. 

Call  this  Problem  4. 

From  (1.7)  we  have 

Rx(A)  = [K  + W(A))/T(A)  = K/T (A)  + 0(A). 

This  suggests  that  it  should  he  possible  to  treat  Problem  4 directly 
and  show  that  a properly  chosen  randomization  among  control  limit  rules 

k 

is  optimal.  For  brevity,  we  will  simply  show  that,  for  properly  etiosen 
values  of  Z,  a control  limit  rule  is  optimal.  Consider  any  fixed  value 
of  K.  We  have  shown  that  among  the  SM  rules  that  minimize  g^  is  (at  least) 
one  of  the  form:  Dispatch  if,  and  only  if,  i^  + i.,  > k.  Let  * be  any  such 

optimal  rule.  If  Z = T(*),  rule  * is  optimal  (among  the  SM  rules)  for 
Problem  4.  To  see  this,  consider  any  SM  rule  A for  which  T(A)  _>  Z.  Then, 
we  obtain  Q(A)  Q(*)  by  noting  that:  g^(A)  _>  g^(*)  or 

Q(A)  + K/T  (A)  _>  Q(*)  + K/T  (*)  or  Q(A)  - Q(*)  _>  K(l/T(*)  - 1/T(A))  >.  0, 
where  the  last  inequality  results  from  T(A)  Z = T(*). 


For  an  example  of  such  analysis  in  an  equipment  replacement  context, 
see  Kolesar  [4]. 
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APPENDIX 


AN  ALTERNATE  APPROACH  TO  THEOREM  I 

We  show  that  an  equivalence  between  Problems  2 and  3 allows  us  to 
claim  that  there  is,  among  the  optimal  SM  rules  for  Problem  2,  one  which 
is  a function  only  of  i^  + i0. 

We  claim  that,  for  any  SM  rule  A for  Problem  2 (our  system) , there 

exists  an  SM  rule  B for  Problem  3,  with  (1.15)  replacing  (1.14)  as  the 
G 

definition  of  p „ , (the  one-dimensional  system)  with  the  following  prop- 
erty: Imagine  an  observer  who  is  informed 

(a)  each  time  an  arrival  occurs  while  the  shuttle  is  waiting 
(but  not  where  it  occurs  for  our  system), 

(b)  each  time  the  shuttle  returns  to  terminal  1,  being  told  the 
total  number  of  passengers  now  waiting,  and 

(c)  each  time  the  shuttle  is  dispatched. 


The  observer  will  not  be  able  to  identify  which  process  is  the  one-dimensional 
system,  controlled  by  B,  and  which  is  ours,  controlled  by  A.  To  see  why, 
consider  the  following  example.  Let  A(0,  0)  = 0,  A(l,  0)  = 1/3,  A(0,  1)  = 4/5, 
and  A(_i)  = 1 for  i^  + 2.  (In  words,  policy  A is:  wait  if  no  one  is 

waiting;  go  if  two  or  more  are  waiting  (anywhere);  go  with  probability  1/3 
if  one  is  waiting  at  terminal  1 and  no  one  at  2;  go  with  probability  4/5 
if  one  is  waiting  at  terminal  2 and  no  one  at  1.)  Then,  the  probability 
that  we  will  dispatch  when  there  is  exactly  one  person  waiting  (somewhere)  is 


B(l) 


= _P0(A1 


r)P0(A2r^ 


~I  ‘ 3 + ~X  ' J + pl(Xlr)pO(X2r/2) 


’0(V)pl(X2r/2)  - l]/[pO(Xlr)pO(X2r/2)  + 

P1(X1r)p0(x2r/2)  + pQ(X;lr)p1(X2r/2)J  . 
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So  the  observer  would  not  be  able  to  tell  our  system  from  the  one-dimensional 
one  that  used  that  value  of  B(l)  (and  B(0)  = 0;  B(k)  = 1,  k >_  2) . 

It  is  easier  to  see  that  there  is  an  equivalent  B for  a given  A than  it 
is  to  specify  it.  Speaking  roughly,  we  do  so  by: 


k 

B(k)  = l 

j=0 


P 


dispatch  from  state  (j,  k - j) 
reach  (j , k - j),  under  rule  A 


reach  state  (j,  k - j) 
under  rule  A 


reach  k passengers' 
waiting  in  total 
under  rule  A 


Formally:  Let 


R(j,  k - j) 


(under  rule  A,  system  enters  \ 
state  (j,  k - j)  before  shuttle! 
is  next  dispatched,  given  that  I 
it  has  just  been  dispatched  / 


and  let 

R(i,  k - i) 

=0 

for  0 ^ j _<  k , k = 0 , 1,  ...  . 

(Explicit  calculation  of  8(j,  k - j)  is  possible  but  sloppy;  it  depends  on  A; 
but  mainly  it  is  irrelevant  to  our  purpose  here.)  Since  a one-step  transition 
to  any  state  is  possible  if^we  dispatch  the  shuttle,  8(j,  k - j)  > 0 for  all 
0 ^ j < k.  By  definition,  f 6(j,  k - j)  = 1 for  k = 0,  1,  ...  . The  equiv- 
alent B for  the  given  A is  ^ ° 

k 

B(k)  = l A( j , k - j)  • 8(1,  k - j)  k = 0,  1 

j=0 


B(j,  k - j)  = R(j,  k - j )/l 


(A . 1) 
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Now,  by  construction  of  the  cost  functions,  if  B is  equivalent  to  A, 
we  have  the  same  average  cost  per  unit  time  for  our  system  and  the  one-dimen- 
sional one.  But,  to  minimize  cost  for  the  one-dimensional  system.  Deb  and 

Serfozo  have  proved  that  among  optimal  rules  (including  non-SM  rules)  is  an 
SM  rule,  call  it  B*,  such  that 

B*(k)  = 0 for  k < k*;  B*(k)  = 1 otherwise.  (A. 2) 

(Actually,  they  treat  a class  of  processes  that  include  (1.3),  (1.4),  and 

(1 . 11)- (1 . 14)  as  a special  case.  An  examination  of  the  proofs  of  their 
results  leading  to  (A. 2)  reveals  that  using  (1.15)  instead  of  (1.14)  does 
not  alter  them.) 

We  claim  that  there  is  an  optimal  SM  rule  for  the  system  (1.1)— (1.5) 
and  (1.9)  which  satisfies  both  (A.l)  and  (A. 2).  Any  rule  A satisfying  (A.l) 
alone  is  equivalent  to  a B which  is  no  better  than  B*,  and  we  now  show  that 
there  is  an  A*  equivalent  to  B*.  It  is  A*(j,  k - j)  = B*(k)  for 

0_<j£k,  k = 0,  1 And,  to  return  to  Problem  1,  we  have  from 

(1.8)  that  an  optimal  rule  for  (1. 1)— (1.5)  and  (1.9)  is  optimal  for 
(1.1)  — (1.6) . That  is,  we  have  proved  Theorem  1. 

OPT  I MAI.  VALUES  OF  THE  CONTROL  LIMIT  FOR  PROBLEM  1 

The  entries  in  the  following  table  are  values  of  k*  such  that,  for  the 
indicated  cost  and  arrival  parameters,  an  optimal  SM  rule  is  to  dispatch  the 
carrier  if,  and  only  if,  a total  of  at  least  k*  passengers  are  waiting  at 
the  two  terminals. 
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Table  A.  1 

OPTIMAL  DISPATCHING  NUMBERS  (k*)  FOR  AN  INFINITE  CAPACITY  SHUTTLE 
OPERATING  WITH  COMPLETE  PASSENGER  WAITING  INFORMATION 

(r  = h = 1) 

Values  of  Values  of  K (trip  cost)  ► 


X1 

0 

1 

2 

4 

8 

16 

32 

X = 

0.25 

0.03125 

1 

1 

i 

1 

2 

3 

4 

0.06250 

1 

1 

i 

1 

o 

3 

4 

0.12500 

1 

1 

i 

1 

2 

3 

4 

0.18750 

1 

1 

i 

1 

2 

3 

4 

0.21875 

1 

1 

i 

1 

2 

3 

4 

u/r 

= 0.25 

0.03125 

1 

1 

i 

2 

3 

4 

5 

0.06250 

1 

1 

i 

2 

3 

4 

5 

0.12500 

1 

1 

i 

2 

2 

3 

5 

0.18750 

1 

1 

i 

2 

2 

3 

4 

0.21875 

- . --  - i 

1 

1 

i 

2 

2 

3 

4 

X = 

0.5 

0.06250 

1 

1 

i 

2 

3 

4 

5 

0.12500 

1 

1 

i 

2 

3 

4 

6 

0.25000 

1 

1 

i 

2 

3 

4 

6 

0.37500 

1 

1 

i 

2 

3 

4 

6 

0.43750 

1 

1 

i 

2 

3 

4 

6 

a/r 

= 0.5 

0.06250 

1 

1 

2 

2 

4 

5 

7 

0.12500 

1 

1 

2 

2 

3 

5 

7 

0.25000 

1 

1 

2 

2 

3 

5 

7 

0.37500 

1 

1 

2 

2 

3 

4 

6 

0.43750 

1 

1 

2 

2 

3 

4 

6 

Values  of 


Values  of 


\ = 4.0  0.5 


1.0 


2.0 


3.0 
3.5 

a/r  =4.0  0,5 

1.0 
2.0 
3.0 


3.5 


X = 8.0  1.0 


2.0 

4.0 

6.0 

7.0 

a/r  = 8.0  1.0 

2.0 

4.0 

6.0 


7 
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Values  of  K (trip  cost) 


-► 


0 

1 

2 

4 

8 

16 

32 

1 

2 

3 

4 

6 

10 

14 

1 

2 

3 

4 

7 

10 

15 

1 

2 

3 

5 

7 

10 

15 

2 

3 

4 

5 

8 

11 

16 

2 

3 

4 

5 

8 

11 

16 

1 

2 

3 

5 

8 

12 

18 

1 

2 

3 

5 

8 

12 

18 

1 

2 

3 

5 

8 

12 

18 

2 

3 

4 

5 

8 

12 

17 

2 

3 

4 

6 

8 

12 

16 

1 

2 

3 

5 

8 

13 

19 

1 

2 

3 

5 

8 

13 

20 

2 

3 

4 

6 

9 

14 

21 

3 

4 

5 

7 

10 

15 

22 

4 

5 

6 

8 

11 

16 

22 

1 

2 

3 

5 

9 

15 

24 

1 

2 

3 

5 

9 

15 

24 

2 

4 

4 

6 

10 

16 

24 

3 

4 

5 

7 

11 

16 

23 

4 


5 


6 


8 


11 
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Values  of 

X1  j 

Values 

0 

; of 
1 

K (trip 
2 

cost ) 
4 

8 

16 

— * 

32 

X = 16.0 

2.0 

1 

3 

3 

5 

9 

16 

25 

4.0 

3 

3 

4 

6 

10 

17 

26 

8.0 

4 

5 

6 

8 

12 

19 

28 

12.0 

7 

7 

8 

10 

14 

21 

30+ 

14.0 

7 

8 

9 

11 

15 

22 

30+ 

ct/r  =16.0 

2.0 

1 

5 

5 

6 

9 

17 

30+ 

4.0 

4 

4 

7 

10 

18 

30+ 

6 

6 

8 

12 

20 

30+ 

8 

8 

10 

14 

21 

30+ 

8 

9 

9 

11 

15 

22 

30+ 
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